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ABSTRACT 

Context. Measuring the water deuterium fractionation in the inner warm regions of low-mass protostars has so far been hampered by 
poor angular resolution obtainable with single-dish ground- and space-based telescopes. Observations of water isotopologues using 
(sub)millimeter wavelength interferometers have the potential to shed light on this matter 

Aims. To measure the water deuterium fractionation in the warm gas of the deeply-embedded protostellar binary IRAS 16293-2422. 
Methods. Observations toward IRAS 16293-2422 of the 53,2 - 44,1 transition of H^**0 at 692.07914 GHz from Atacama Large 
Millimeter/submillimeter Array (ALMA) as well as the 3i.3 - 22.0 of H^**0 at 203.40752 GHz and the 3i.2 - 22.1 transition of HDO at 
225.89672 GHz from the Submillimeter Array (SMA) are presented. 

Results. The 692 GHz H^^O line is seen toward both components of the binary protostar. Toward one of the components, "source B", 
the line is seen in absorption toward the continuum, slightly red-shifted from the systemic velocity, whereas emission is seen off-source 
at the systemic velocity. Toward the other component, "source A", the two HDO and H^^O lines are detected as well with the SMA. 
From the H'^^O transitions the excitation temperature is estimated at 124 ± 12 K. The calculated HDO/H2O ratio is (9.2 ± 2.6) x 10"'* 
- significantly lower than previous estimates in the warm gas close to the source. It is also lower by a factor of ~5 than the ratio 
deduced in the outer envelope. 

Conclusions. Our observations reveal the physical and chemical structure of water vapor close to the protostars on solar-system 
scales. The red-shifted absorption detected toward source B is indicative of infall. The excitation temperature is consistent with the 
picture of water ice evaporation close to the protostar. The low HDO/H2O ratio deduced here suggests that the differences between 
the inner regions of the protostars and the Earth's oceans and comets are smaller than previously thought. 

Key words, astrochemistry - stars: formation - protoplanetary disks - ISIVI: abundances - ISIVI: general 



1. Introduction 
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Water plays an essential role for life as we know it, but its ori- 
gin on Earth is still unclear: was water accreted during the early 
stages of Earth's formation, or brought by smaller solar system 
bodies such as comets at later times? To deduce the origin of 
Earth's water and the amount of chemical processing it has ex- 
perienced, one option is to measure the water deuterium frac- 
tionation (HDO/H2O) during different stages in the evolution of 
protostars and compare it to what we measure in Earth's oceans 
and comets. 

G enerally the HDO/Hp O ratio in Earth's oceans of 3 x 10 
(e.g., Lecuyer et al. 1998 1 and Oort cloud comets of 8.2x10 
( Villan ueva et aL 2009| l are found to be enhanced above the Dfl^ 
ratio in the protosolar nebula ~ 1 .5 x 1 0"*" ( Linsky 2003 ; Geiss & 
|Gloeckler||1998| l due to deuterium fractionation processes. The 
factor of 2 higher abundance ratio in the Oort cloud comets than 
in Earth's oceans has previously been taken as an indication that 
only a small fraction of Earth's water was delivered by comets. 
Recently however, a HDO/H2O ratio of 3.2 x 10 was measured 



A HDO/H2O ratio of "x" corresponds to a D/H ratio of "x/2". 



for the Jupite r class comet Hartle y 2 with the Herschel Space 
Observatory ( [Hartogh et al.||201l| and 4.2 x 10 "^ f or the Oort 
cloud comet Garradd ( [Bockelee-Morvan et al.|2012| i, indicating 
values closer to those of Earth's water. 

Attempts at measuiing the water deuteiium fractionation in 
protostars have resulted in different conclusions. Parise et al. 
(|2003 ) used ground-based infrared observations of the stretch- 
ing bands of OH and OD in water ice in the outer parts of en- 
velopes and found upper limits ranging from 0.5% to 2% for the 
HDO/H2O ratios in four embedded low-mass protostars. In the 
gas-phase it is possible to detect lines of HDO, but such studies 
differ on the interpretation with HDO/H2O ratios in protostars 
ranging from cometary values (Stark et al. 2004 j ), to a few % 
( Parise et al.|2005t|Liu et al.|201 l| l. Even more recently, [Coutens] 
let al.| ( |2012 !) deduced a HDO/H2O ratio in IRAS 16293-2422 of 
3.4 X lO^in the inner parts and 0.5 x 10"^ in the outer envelope 
by modeling a large range of lines observed with Herschel. 

One problem with previous measurements of HDO/H2O is 
the relatively large beam size of single-dish ground- and space- 
based telescopes. Spheiically symmetric power-law models of 
protostellar envelopes have usually been employed to interpret 
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the observations. While such models are appropriate to interpret 
continuum and line emission on larger scales (>300 AU), they 
are not suited to unambiguously analyze the observed compact 
components since there are clear indications that they are not an 
accurate representation of the conditions on small scales (e.g., 
J0rge nsen et al.|[2005 Chiang et al.||2012] l. Estimates of abun- 
dance ratios on these smaller scales are thus subject to signifi- 
cant uncertainties due to extrapolations of the underlying physi- 
cal structures. 

High angular resolution millimeter wavelength aperture syn- 
thesis observations offer a possibility to circumvent this issue. 
Recently J0rg ensen & van Dishoeck| |201Q) detected the water 
isotopologue Hj^O toward the deeply embedded protostar NGC- 
1333 1RAS4B on scales of <50 AU using the IRAM Plateau de 
Bure Interferometer (PdBl), which combined with an upper limit 
on the HDO column density from the SMA resulted in a 3cr up- 
per limit to the HDO/H2O abundance ratio of 6x 10 To follow- 
up these results we initiated an extended survey of the Hj^O and 
HDO emission on arcsecond scales using the IRAM PdBI and 
SMA ( iPerTson et al.|2012l l. 

IRAS 16293-2422 is a Class protostellar binai-y (sep ~ 5", 
600 AU) located 120 pc away in the LDN 1689N cloud in the 



p Ophiucus star-forming region (Knude & H0g 



1998 Loinard 



|et al.|[2008]). With a rich spectrum at (sub)millimeter wave 
length s (Blake et al.||19 94; van Dishoeck et al. i'1995; Cazaux' 
et al.||2003, Chandler et al.,2005j " ]Caux et al.||2011, J0rgensen| 
et al.(|20l"l| ) it has been one of the prime targets for studies 
of astrochemistry during the star-formation process, revealing 
the presence of a range of complex organic species (BottinelH] 
^taLf 2004 , Ku an et al.]2 004 ; Bisschop et al. 2008 ) and prebiotic 
molecules (J0rgensen et al.^2012) on (sub)arcsecond scales. 

In this letter we present, for the first time, high-resolution 
ground-based observations of several isotopologues of water to- 
ward IRAS 16293-2422 with both the ALMA (~0'.'2; 24 AU) 
and the SMA (~2'.'3; 276 AU). We derive direct, model- 
independent estimates of the water excitation temperature and 
use it to calculate the column density and the HDO/H2O ratio in 
the warm inner envelope. 



2. Observations 

Observations of the deeply-embedded low-mass protostellar bi- 
nary IRAS 16293-2422 were carried out at 690 GHz with 
ALMA and 230 GHz with the SMA, tai-geting the 53,2 - 44j 
(692.07914 GHz) and 3i 3 - 220 (203.40752 GHz) transitions of 
HfO and the 3i 2 -22 1 (225.89672 GHz) transition of HDO (see 



Table A.l in the appendix/electronic version). In the tables, the 
203 GHz observations are indicated with "1", 225 GHz observa- 
tions with "2" and 692 GHz with "3" (see Tableland TablejAl) 
The ALMA observations of IRAS 16293-2422 were con- 
ducted as part of the ALMA Science Verification (SV) pro- 
gram: IRAS 16293-2422 was observed with 13 antennas on 
April 16 and 17, 2012 in a seven pointing mosaic centered at 
a = 16''32'"2Z7, S = -24°28'32'.'5 (J2000). The observations 
cover projected baselines from 26 to 402 m (62 to 945 k/l). One 
of the basebands was centered at 691.299 GHz with a band- 
width of 1.875 GHz and a spectral resolution of 0.923 MHz 
(0.4 km s"'), a setup that covers the Hj^O 53,2 - ^a,\ line at 
692.07914 GHz. Calibration observations include the quasars 
1924-292 and 3c279 for the bandpass, the asteroid Juno for 
the amplitude, and the quasars 1625-254 and nrao530 for the 
phase. The science verification data are available as calibrated 
Mv-data sets, which were used in our analysis. The calibrated 



data were imaged using the CASA software package ( |McMullm| 
let al.|2007| l. 

The lower-lying excited Hj^O 3 1,3 - 22,0 transition at 
203.4075 GHz was observed with the SMA on May 1, 2011 in 
the compact configuration with seven antennas. This configura- 
tion resulted in projected baselines between 9 and 69 meters (6 
to 47 k/l). The passband was calibrated by observations of the 
quasar 3c279 while the absolute flux and complex gains were 
calibrated by observing Titan and the quasars 3c279, 1517-243 
and 1626-298. The raw data calibration followed the standard 
r ecipes using the MIR package (Qi 2008) and then MIRIAD 
(Saul t et al.|199 5) was used to subtract the continuum from the 
data to create continuum-free line maps. 

Finally, we utilized SMA observations of the HDO 3i_2 - 
22,1 line at 225.89672 GHz from the SMA ('J0rgensen et al 



|201 l| l: those observations have a spectral resolution of 0.41 MHz 
(0.54 km s"') and cover the projected baselines between 8.6 - 
1 19.6 m (6.5-90 k/i). For further information about those obser- 
vations we refer to |j0rgensen et al. ( |201 l| l. The resulting beam 
size, field of view, velocity resolution and RMS of the various 



line data are summarized in Table A. 1 in the appendix/electronic 



3. Results 

Figure [T] shows the spectra around the water lines toward the 
continuum peaks with Gaussian fits. Source A is in itself a bi- 
nary ( [Chandler et al. 2005) and we here refer to the components 
as Al (the Northeast component) and A2 (the Southwest com- 
ponent) and the lower resolution single component as A. Both 
Hj^O lines are clearly detected toward source A in emission, 
while only the 692 GHz Hj^O line is seen toward source B in 
absorption. The H2^0 emission is clearly associated with peak 
Al in the 692 GHz data. The H^'^O 692 GHz absorption line to- 
ward the continuum peak of source B is narrow and marginally 
red-shifted; the absorption is not due to the broad outflow seen 
in some lines (I J0rgensen et al.|201 l| l. 

This weak absorption feature is not seen in the SMA data 
but the high sensitivity of ALMA combined with the stronger 
continuum at high frequencies is enough to detect it. The lines 
were identified using the JPL and CDMS catalogs through the 
Splatalogue compilation ( |Pickett et al.|1998tpUller et al.|20"0T] l. 

The Splatalogue compilation shows no other possible line 
candidates at least +1 km s"' from the line peak positions of the 
H^'^O line at 203 GHz and the HDO hne at 225 GHz. Lines that 
fall within ±2-3 km s"^ exist, but are of too low intrinsic line 
strengths or lack additional components that should have been 
detected at other velocity offsets e.g., (CH3)2CO, '''CH3CH2CN 
or CgH. At 203.403 and 203.410 GHz (i.e., about 10 and - 
3 km s"') lies CH3OCH3 33,1 - 22.1EA/EE, which is accounted 
for when fitting the 203 GHz line. These lines are narrow enough 
to avoid any significant interference with Hj^O (see top most plot 
inFig.[T]i. 

For the H2^0 ALMA data at 692 GHz the line identification 
in source A is complicated by the fact that the source is resolved 
and shows a systematic velocity pattern in other lines ([Pineda] 
et al. 2012). The line toward the continuum position of source 
A is red-shifted by about 2 km s"' from the systemic i^i; » of 



3.2 km s~', similar to lines from other species (see Fig.fTl panel 
four). The third spectrum from the top in Fig. TJshows the spec- 
trum after smoothing with a Gaussian kernel of 1" size. The 
Gaussian fit to these data peaks closer to visr-^-'^ km s"', in- 
dicative indeed of resolved velocity structure. Outflow emission 
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(c.f., |Loinard et al.|2012[ ), as also seen in H2O masers on small 
scales ( Alves et al.|2012| ), could explain this velocity shift, but it 



is not possible to deduce its importance from the current data 

However, another possibility is that line blending causes the 
shift. Toward source B, the off position spectrum shows two 
emission lines, the second line at ~ 6 km s ' is unidentified (bot- 
tom spectrum Figjlb. The H^^O line at 692 GHz toward source 
A may therefore also be affected by line blending. However, the 
unambiguous detection of the 692 GHz line toward source B and 
the SMA observations of water toward source A suggest that at 
least half of the 692 GHz feature toward source A is indeed due 
to the H^^^O line. 

Table [1] lists parameters from Gaussian fits to the image and 
uv-plane of the line emission. The integrated intensity maps for 
the detected water lines are shown in Fig. [2] 
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Fig. 1. Spectra of the targeted water lines toward both IRAS 16293- 
2422 A and B. The 203 GHz Hj'^O line spectrum (top) was binned to 
twice the resolution for clarity. The ID Gaussian fits (red: H^^^O; green: 
CH3OCH3; magenta: Unidentified) and the RMS (blue) are plotted. 
Shaded areas show the interval over which the integrated intensities 
have been calculated (yellow and blue). The dotted green vertical lines 
shows ulsr = 3.2 and 2.7 km s"' for source A and B respectively. The 
third spectrum from the top (blue fill) is from data smoothed with a 1" 
Gaussian kernel. 



4. Discussion 

4.1. Spectra and integrated intensity 

The different water lines toward source A show similar charac- 
teristics. The peak positions and line widths of the Gaussian 
fits to the 203 GHz H'^O and 225 GHz HDO lines are simi- 
lar and agree with previous studies of other molecules toward 
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Fig. 2. Integrated intensity maps for all observed water lines calculated 
from channels ±3 from the systemic velocity, deduced from previous 
observations ( J0rgensen et al. 201 1 1. The absorption line toward source 
B was integrated in the interval 0-5 km s"'. Note the different spatial 
scales of the top vs the bottom panels. The beam is shown in the lower 
right corner and the gray crosses show the position of continuum peaks 
from elliptical Gaussian fits. Units on the axes are offset in arcseconds 
from the phase center of the 203 GHz observations. Contours are in 
steps of 2cr starting at 3cr, dashed contours represent negative values. 

Table 1. Parameters from fits to the integrated maps and spectra of the 
different water lines for source A and B. The given errors are the statis- 
tical uncertainties and the 20% calibration uncertainty in intensity is not 
included. The spectral line fits were done toward the continuum peak 
and integrated intensities deduced from circular Gaussian fits in the uv- 
plane (203 and 225 GHz) or an elliptical Gaussian fit in image plane 
(692 GHz). 



Line 


Size 


Intensity 


Line width 


Id 


["] 


[Jy km s"'] 


[kms-i] 


1(A) 


1.7 ±0.7 


3.9 ± 1.4 


6.9 ± 1.5 


2(A) 


1.1 ±0.4 


8.9+1.8 


7.0 ± 0.7 


3(A1) 


0.2x0.4(13°) 


3.4 ± 0.6 


6.3 ±0.7 


3(B) 


Point fit 


-0.43 ± 0.07 


0.9 ± 0.2 



Notes. The column "Line Id" gives the source indicated in parentheses, 
and the number corresponds to one of the lines (see text and Table |A. l| . 
The line width is the FWHM from a Gaussian fit to the spectral line. 
The elliptical Gaussian parameters (for line 3 (Al)) are given as minor 
and major axis (±0'.'l) and position angle (PA, ±4°). 



IRAS 1 6293-2422 e.g.,|Bisschop et aL]([2OO8ll; |j0rgensen et al 
( |20TT] l; |j0rgensen et al.| ( |2012| l; [Pineda et aL] ( |2012| ) (seTTV 
ble yj. The mapped emission is compact and traces the warm 
water on scales 7?<200 AU. 

Toward source B, the 692 GHz Hj^O line is detected in ab- 
sorption. The spectral line is slightly red-shifted compared with 
=2.7 km s"' of source B, consistent with the picture of on- 



vlsr- 



going infall in this source (j0rgensenetaL 2012; Pineda et al. 
|2012) l. Integrating over a larger velocity interval reveals, in addi- 
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tion to the absorption, a blue-shifted emission peak, offset ~ 0'.'2 
south of the absorption (see Fig.[T]and[2]). 



4.2. Temperature and deuteration of water vapor 

To calculate the excitation temperature and column densities to- 
ward source A we assume LTE conditions and that the emis- 
sion at 203 and 225 GHz has the same extent as in the 692 GHz 
observations (see Table [IJ. Scaling the 203 GHz Hj**© intensi- 
ties to the 692 GHz source size gives an excitation temperature 
of Tex = 124 ± 12 K for Hf O. If the 692 GHz Hne is indeed 
blended, it would imply that the intensity estimate, and also the 
excitation temperature, are upper limits. Halving the intensity 
of the 692 GHz line causes a drop to T^^ = 107 + 9 K. The ex- 
tent and peak location of the emission in source A show that the 
warm water vapor is located at a small projected distance from 
the center of the protostar. The determined excitation temper- 
ature indicates that the distance is also small along the line of 
sight. 

The estimated HDO column density is 4.8 x 10'^ cm"^, cor- 
rected for beam dilution when assuming that the extent of the 
emission is that of the 692 GHz Hj^O observations and Tex - 
124 K. The gas-phase H2O column density is 5.3 x 10^" cm"^ 
assuming the same as for HDO and that the isotopic abundance 
ratio of '*0/'**0 is 560 in the local interstellar medium (Wilson, 
|& Rood|1994| l. Assuming an uncertainty of about 20% for the 
column densities, originating in the flux calibration, the best es- 
timate of the HDO/H2O ratio is (9.2 + 2.6) x 10 1 

Given the uncertainty in the determination of the excitation 
temperature, testing the effect of different temperatures is im- 
portant. If the excitation temperature is as low as 80 K the 
HDO/H2O ratio becomes 7.8 + 2.2 x 10 Increasing the ex- 
citation temperature to 300 K P0rgensen et al.|2012) i increases 
the ratio to 1.1 ± 0.3 x 10"^. The conclusions in this paper do 
not change over this wide interval in Tex. In addition, since the 
203 GHz H^^O and 225 GHz HDO observations have compa- 
rable beam and sources size, and arise from levels with similar 
energies, the deduced HDO/H2O ratio is robust. However, if the 
HDO emission were more extended, its column density and the 
HDO/H2O ratio would be lower In this scenario our inferred 
ratio is an upper limit. 

In contrast with previous estimates based on models of 
single-dish observations (Par ise et al.|20 05 ; Cout ens et al.|2 012 ), 
the deduced HDO/H2O ratio in flie warm gas of IRAS 16293- 
2422 is only slightly higher than found in Earth's oceans and by 
recent Herschel observations of comets. Given the possible sys- 
tematic errors due to assumptions of the extent, they could be 
even closer. Within the statistical uncertainties our observed ra- 
tio for this protostar agrees with the earlier ratios for Oort cloud 
comets. Comparing these different ratios for water directly as- 
sumes that the reservoirs, i.e., comets. Earth (planet) and inner 
protostellar region are linked and can be related. 

The difference with the earlier estimates comes from the fact 
that those data are sensitive to much larger scales than the high- 
resolution interferometric observations presented here, which di- 
rectly image the water emission on 25 - 280 AU scales. Our in- 
terferometric observations provide a strong, model independent 
constraint on the deuteration of water in the innermost regions 
of protostars. That the low HDO/H2O ratio in the warm gas is 
not much different from the cometary values is an indication that 
significant processing of the water between these early stages 
and the emerging solar system is not required. 



Further high-resolution interferometric measurements to- 
ward larger samples of protostars will reveal whether the warm 
HDO/H2O ratio is similar in different protostars. In particular, 
future high angular resolution observations with ALMA will be 
able to resolve possible variations in the HDO/H2O ratio with 
distance from the central protostar and thereby show whether the 
slightly different ratios measured in different types of comets po- 
tentially could be related to their spatial origin in the protostellar 
envelope ( Robert et al.|2000 1. 
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Table A.l. Relevant parameters from the molecular line catalogs for the observed lines. 



Line 


Species 


Frequency 


Transition 


Line strength 




Beam 


Resolution 


RMS 


Id' 




[GHz] 




[Debye^] 


[K] 


Size (PA) 


[kms-i] 


['] 


1 




203.40752 


3 1,3 - 22,0 


0.34 


203.7 


3.1" X 5.0" (-25.3°) 


0.30 


208 


2 


HDO 


225.89672 


3l,2 - 22,1 


0.69 


167.6 


2.3" X 3.9" (15.0°) 


0.54 


190 


3 


H>«0 


692.07914 


53,2 - 44,1 


1.26 


727.6 


0'.'17x0'.'29 (109.8°) 


0.40 


82 



' Identifies the lines in Table [T] 
mJybeam"' channel ' 



